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Abstract — This paper studies the achievable rate region of 
the if -user Gaussian multiple-input single-output interference 
channel (MISO-IC) with the interference treated as noise, when 
improper or circularly asymmetric complex Gaussian signaling 
is applied. The transmit optimization with improper Gaussian 
signaling involves not only the signal covariance matrix as in the 
conventional proper or circularly symmetric Gaussian signaling, 
but also the signal pseudo-covariance matrix, which is convention- 
ally set to zero in proper Gaussian signaling. By exploiting the 
separable rate expression with improper Gaussian signaling, we 
propose a separate transmit covariance and pseudo-covariance 
optimization algorithm, which is guaranteed to improve the 
users' rates over the conventional proper Gaussian signaling. In 
particular, for the pseudo-covariance optimization, we establish 
the optimality of rank-1 pseudo-covariance matrices, given rank- 
1 transmit covariance matrices for achieving the Pareto boundary 
of the rate region. Based on this result, we significantly reduce 
the number of variables in the pseudo-covariance optimization 
and thereby develop an efficient algorithm for it by applying 
the celebrated semidefinite relaxation (SDR) technique. Finally, 
we extend the result to the Gaussian MISO broadcast channel 
(MISO-BC) with linear precoding at the transmitter. 

Index Terms — Improper Gaussian signaling, interference chan- 
nel, beamforming, pseudo-covariance optimization, semidefinite 
relaxation, broadcast channel, linear precoding. 



I. Introduction 

The interference channel (IC) models the practical scenario 
in wireless communication when more than one transmitter- 
receiver pairs communicate independent messages at the same 
time and thus interfere with each other due to the shar- 
ing of a common spectrum. Characterizing the information- 
theoretic limit for the general if -user IC is a long-standing 
open problem | |1|2[ . For the two-user single-input single- 
output IC (SISO-IC) with additive Gaussian noise, it was 
recently proved that one particular form of the celebrated 
Han-Kobayshi scheme achieves within one bit to the capacity 
region |3l. Since the capacity-approaching scheme in general 
requires multi-user encoding and decoding, which are difficult 
to implement in practical systems, a great deal of research 
on Gaussian ICs has focused on characterizing the achievable 
rate regions ||4)-|fT2"), under the assumption of employing the 
single-user decoding (SUD) with the interference treated as 
Gaussian noise at receivers. 
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For example, the achievable rate region with the interference 
treated as noise has been completely characterized for the 
Gaussian SISO-IC Q, single-input multiple-output IC (SIMO- 
IC) (5), and multiple-input single-output IC (MISO-IC) (6)- 
(9). For the more general Gaussian multiple-input multiple- 
output IC (MIMO-IC), some recent results were also reported 
in ||l0|-p2). In general, the achievable rate region of an IC 
is completely characterized by its Pareto boundary, which 
constitutes all the achievable rate-tuples for all users at each 
of which it is impossible to improve one user's rate, without 
simultaneously decreasing the rate of at least one of the other 
users. A traditional approach for such a characterization is 
via solving a sequence of weighted-sum-rate maximization 
(WSRMax) problems AS). An alternative method based on the 
concept of "rate profile" was proposed in |7J for the MISO-IC 
setup, which eventually results in solving a sequence of signal- 
to-interference-plus-noise ratio (SINR) feasibility problems 
that are easier to handle than WSRMax problems. It is worth 
mentioning that for the MISO-IC setup, it has been shown 
by various methods in [7|-|9| that all the rate-tuples on the 
Pareto boundary can be achieved with transmit beamforming, 
i.e., with rank-1 transmit covariance matrices. 

Besides the achievable rate region characterization, signif- 
icant research effort on Gaussian ICs has been devoted to 
maximizing the (weighted) sum-rate achievable by all the users 
| 5|13| -p6|. Unfortunately, such problems have been shown to 
be NP-hard in general 1 15|16| . Many suboptimal or heuristic 
algorithms have thus been proposed, e.g., the gradient descent 
algorithm [17], the interference-pricing based algorithm |18|, 
the game-theory based algorithm |19| , the egoistic versus 
altruistic based algorithm ]20[ , and the iterative weighted min- 
imum mean-square-error (MMSE) based algorithm |21|22] . 
More recently, for Gaussian SISO-IC, SIMO-IC and MISO- 
IC, the globally optimal solutions to WSRMax problems have 
been obtained under the monotonic optimization framework 
1 5 23 1-[26 1. However, the complexity of such globally optimal 
algorithms increases exponentially with the number of users, 
and their generalization to the more general MIMO-IC remains 
unknown. Moreover, it is worth mentioning that there have 
been a great deal of research interests over the last few years 
in studying Gaussian ICs from the degrees-of-freedom (DoF) 
perspective by using the technique known as interference 



alignment (IA) (see [27| and references therein). 



However, the aforementioned works on ICs have mostly 
assumed proper or circularly symmetric complex Gaussian sig- 
naling for transmitted signals. It was first revealed in |28| that 
the more general improper or circularly asymmetric complex 
Gaussian signaling, jointly applied with symbol extension and 
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IA, is able to improve the achievable DoF for the sum-rate of a 
three-user Gaussian SISO-IC at the asymptotically high signal- 
to-noise ratio (SNR). Later, it was shown that even for the two- 
user SISO-IC where IA is not applicable, the achievable rate 
region can still be enlarged with improper Gaussian signaling 
over the conventional proper Gaussian signaling [ 29|30| . Re- 
cently, we have shown in our previous work ]3l| that with 
improper Gaussian signaling, the user's achievable rate in the 
general MIMO-IC can be expressed as the summation of the 
rate achievable by the conventional proper Gaussian signaling, 
which depends on the users' transmit covariance matrices only, 
and an additional term, which is a function of both the users' 
transmit covariance and pseudo-covariance matrices. Such a 
separable rate structure was exploited in pT) to optimize 
the covariance and pseudo-covariance separately so that the 
obtained improper Gaussian signaling strictly outperforms 
the conventional proper Gaussian signaling in terms of the 
achievable rate region. However, the algorithms proposed in 
pi) are for the two-user SISO-IC and cannot be applied when 
there are more than two users and/or multiple transmitting 
antennas. This thus motivates our current work that extends 
the result in J3T) to optimize improper Gaussian signaling to 
the more general K-usei MISO-IC. 

Similar to J7J, we apply the rate-profile technique to char- 
acterize the achievable rate region of the MISO-IC with the 
interference treated as Gaussian noise. However, unlike the 
case with proper Gaussian signaling, the resulting optimization 
problem with improper Gaussian signaling is shown to be non- 
convex and thus difficult to be solved optimally. By adopt- 
ing the similar separate covariance and pseudo-covariance 
optimization approach as in [31] , we develop an efficient 
suboptimal solution for this problem in the MISO-IC setup. 
Specifically, for the pseudo-covariance optimization, we first 
establish the optimality for rank- 1 pseudo-covariance matrices, 
given rank-1 transmit covariance matrices for achieving the 
Pare to optimal rates. Moreover, we show that each rank- 
1 pseudo-covariance matrix is parameterized by one sin- 
gle complex scalar. Based on this result, we formulate the 
original matrix optimization problem to an equivalent vector 
optimization problem in much lower dimensions. We then 
apply the celebrated semidefinite relaxation (SDR) technique 
1 32 33 1 to find an efficient approximated solution for the 
reformulated problem. It is worth noting that the approach of 
using SDR for solving non-convex quadratically constrained 
quadratic programs (QCQPs) has been successfully applied to 
find high-quality approximated solutions for various problems 
in communications and signal processing (see [32 33] and 
references therein). For our pseudo-covariance optimization 
problem in the if -user MISO-IC setup, we show that the 
proposed SDR-based solution is indeed optimal when K = 2. 
Finally, we show that the improper Gaussian signaling scheme 
developed for the MISO-IC can also be applied to the if -user 
MISO broadcast channel (MISO-BC), under the assumption of 
employing linear precoding at the multi-antenna transmitter. 

The rest of this paper is organized as follows. Sec- 
tion [n] presents the system model and problem formulation. 
Section III develops the separate covariance and pseudo- 



proposed algorithm is extended to the MISO-BC with linear 
precoding. Section [V] presents numerical results. Finally, we 



conclude the paper in Section VI 



Notations: In this paper, scalars are denoted by italic letters. 
Boldface lower- and upper-case letters denote vectors and 
matrices, respectively. For a square matrix S, Tr(S) denotes 
the trace. S > and S >- mean that S is positive 
semidefinite and positive definite, respectively. C MxN and 
jjMxiv (jgjjotg tjjg S p ace of MxN complex and real matrices, 
respectively. For an arbitrary matrix A, A*, A T , A H and 
rank(A) represent the complex-conjugate, transpose, Hermi- 
tian transpose and rank of A, respectively. The symbol i 
denotes the imaginary unit, i.e., i 2 = — 1. [v]j denotes the jth 
element of the vector v, while 1 1 v | denotes its Euclidean norm. 
For a complex number x, \x\ denotes its magnitude. CAf(x, S) 
represents the circularly symmetric complex Gaussian (CSCG) 
random vector (RV) with mean x and covariance matrix S. 
5ft{-} and represent the real and imaginary parts of a 

complex number, respectively. 

II. System Model and Problem Formulation 

This paper considers a if-user MISO-IC, where each trans- 
mitter is intended to send independent information to its 
corresponding receiver, while possibly interfering with all 
other K— 1 receivers. Assume that each transmitter is equipped 
with M > 1 antennas and each receiver with one single 
antenna. The received baseband discrete-time signal for user 
k is given by 

y k (n) = h kk x k (n) + h kj Xj(n) + v k (n), Vfc, (1) 



where n is the symbol index; h kk £ C lx denotes the 
direct channel from transmitter k to receiver k, while h k j E 
C lxM ,j k, denotes the interference channel from trans- 
mitter j to receiver k; we assume quasi-static fading and 
thus all channels are constant over n's for the problem of 
our interest; v k {n) represents the independent and identically 
distributed (i.i.d.) zero-mean CSCG noise with variance a 2 , 
i.e., v k (n) ~ CJ\f(0, a 2 ); and x k (n) € C Mxl is the transmit- 
ted signal vector from transmitter k, which is independent of 
Xj (n) for j 7^ k. In this paper, for the purpose of exposition, 
we assume that the technique of symbol extensions over 
time [28] is not used. Hence, xj,(n) is independent over n. 
For brevity, n is omitted in the rest of this paper. Different 
from the conventional setup where proper Gaussian signaling 
is assumed, i.e., x^ ~ CAf(0, C Xfc ), with C Xfc denoting 
the transmit covariance matrix, in this paper, we consider 
the more general improper Gaussian transmitted signals. For 
the background knowledge of improper (Gaussian) random 
vectors (RVs), the readers may refer to |31|34| and references 
therein. 

For the zero-mean transmitted Gaussian RV Xfe, we denote 
its covariance and pseudo-covariance matrices as C Xfc and 
C Xfc , respectively, i.e., 



C Xfc =E(x fe xf), C Xfc =E(x fe x£), 



(2) 



covariance optimization algorithm in details. In Section IV the 



where C Xfc is Hermitian and positive semidefinite, and C Xfc 
is symmetric in general. For the conventional proper Gaussian 
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signaling, the pseudo-covariance matrices C Xfc 's for all trans- 
mitters are set to zero matrices, and thus are not included in the 
transmit optimization. However, for the more general improper 
Gaussian signaling considered in this paper, the additional 
degrees of freedom given by the pseudo-covariance matrices 
provide a further opportunity for improving rate over proper 
Gaussian signaling |27|. C Xfc and C Xfc are a valid pair of 
covariance and pseudo-covariance matrices, i.e., there exists a 
RV Xfc with covariance and pseudo-covariance matrices given 
by C Xfc and C Xfc , respectively, if and only if the corresponding 
augmented covariance matrix C Xfc defined below is positive 
semidefinite |35] , 



r A 



y o. 



(3) 



For the MISO-IC with single-antenna receivers, the covariance 
and pseudo-covariance of the received signal y k can be written 

as 



K 



C, fe =E(j /fe2/ *)=^h fcj C Xj hf 3+( T 2 , 



K 



c y k = K (ykVk) = 5Z hfe J Cx i h 



(4) 



(5) 



Denote the interference-plus-noise term at receiver k as s&, 
i.e., s k = Yljatk h/y Xj + v k . Then the covariance and pseudo- 
covariance of s k are given by 

C Sk = E(s kS * k ) = ]T ^ S 'x K + °\ (6) 
C Sk =E( SkSk ) = J2^C X] hl r (7) 



Under the assumptions of Gaussian inputs and that the in- 
terference is treated as improper Gaussian noise at receivers, 
and by applying the result in pT) to the MISO-IC setup, the 
achievable rate at receiver k can be expressed as 



Rk - log 1 



hA; fe C Xfc hf fc 



° 2 + E^ k hfciCx,hg 



-log- 

2 6 i-c s - 2 iaJ 2 



(8) 



It is observed from ^ that with improper Gaussian signal- 
ing, each user's achievable rate is a summation of the rate 
achievable by the conventional proper Gaussian signaling, 
denoted by i?^. loper ({C x }), and an additional term, which 
is a function of both the users' covariance and pseudo- 
covariance matrices. Therefore, for a given set of transmit 
covariance matrices obtained by any proper Gaussian signaling 
scheme, the achievable rates in MISO-IC can be improved 
with improper Gaussian signaling by choosing the pseudo- 
covariance matrices that make the second term in |8]l strictly 
positive. 

The achievable rate region 7Z for the K-nsei MISO-IC is 
defined as the set of rate-tuples that can be simultaneously 



achieved by all users under a given set of transmit power 
constraints at each transmitter, denoted by P k , k = 1, • • • , K. 
With R k given in we thus have 



u 



{n,--- ,r K )-0<r k <R k ,Vk\, (9) 



Tr(C Xfc )<P fc , 
C x ^0, V/c 



where the constraint C Xfc y follows from (j3). 

To characterize the Pareto boundary of the achievable rate 
region 7Z, we adopt the rate-profile method as in J7]. Specif- 
ically, any Pareto-optimal rate-tuple on the boundary can be 
obtained by solving the following optimization problem with 
a given rate-profile vector denoted by a— (a\ ■ ■ ■ ax)- 

(PI): max. R 

{c x , },{c x . 



s.t. R k > a k R, \/k, 

Tr(C x J <P k , Vfc, 

y o, Vfc, 



where a k denotes the target ratio between user fc's achievable 
rate and the users' sum-rate, R. Without loss of generality, we 
assume a k > 0,Vfc, and Ylh=x ak = Denote the optimal 
value of (PI) as R*. Then the rate-tuple i?*- a must be on 
the Pareto boundary of the rate region 7Z. Thereby, by solving 
(PI) with different rate-profile vectors a, the complete Pareto 
boundaiy of 1Z can be found. 



III. Separate Covariance and Pseudo-Covariance 
Optimization 

(PI) is a non-convex optimization problem, and thus it 
is difficult to achieve the global optimum efficiently. In 
this section, we propose a separate covariance and pseudo- 
covariance optimization algorithm by utilizing the separable 
rate expression given in ([8]) to obtain an efficient suboptimal 
solution for (PI). Specifically, the covariance matrices of the 
transmitted signals are first optimized by setting the pseudo- 
covariance matrices to zero, i.e., proper Gaussian signaling is 
applied. Then, the pseudo-covariance matrices are optimized 
with the covariance matrices fixed as the previously optimized 
values. With such a separation approach, the obtained im- 
proper Gaussian signaling design is guaranteed to improve 
the achievable rates over the conventional proper Gaussian 
signaling. 



A. Covariance Optimization 

When restricted to proper Gaussian signaling by setting 
C Xfc = 0,V/c, the second term in ([8]) vanishes to zero and 

If a^ Q = for any ko £ {1. • • • , K} , then we can simply ignore user fco 
in the optimization problem (PI) by removing its corresponding constraints. 
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(PI) reduces to 
(Pl.l): max. r 



s.t. log 1 



h fefc C Xfc 1 1 /. /. 



> a k r, Vfc, 



■ y. J ,iM,c X! h>> l 

Tr(C x J <P k , Vfc, 
C Xfe t 0, Vfc. 

Denote the optimal value of (Pl.l) as r* . Then the rate-tuple 
r* ■ a is on the Pareto boundary of the achievable rate region 
with proper Gaussian signaling. It has been shown in J7|-||9} 
that all the Pareto-optimal rate-tuples with proper Gaussian 
signaling can be achieved by rank-1 transmit covariance ma- 
trices or beamforming. Therefore, without loss of optimality 
for (Pl.l), let 



C Xfc =t fc tf, Vfc, 



(10) 



where t k is the transmit beamforming vector for user fc. Then 
for any fixed target rate r, the feasibility problem related to 
(Pl.l) can be formulated as 



(P1.2): Find {t fc } 

s.t. a 2 + J2 Hjtj\ 2 < — _ ^ hfcfctfc) 2 



Vfc, 



%{h kk t k } = 0, \\t k \\ 2 <P k , Vfc, 

where without loss of generality, we have assumed that for 
each user fc, h. kk t k is a nonnegative real number J36) . (PI. 2) 
is a second-order cone programming (SOCP) problem, which 
can be efficiently solved (37j. If (PI. 2) is feasible, then the 
optimal value of (Pl.l) satisfies r* > r; otherwise, r* < r. 
Therefore, (Pl.l) can be optimally solved by solving (PI. 2) 
with different values of r, and applying a bisection search over 
r (371. 



B. Pseudo-Covariance Optimization 

Denote the optimal solution to (Pl.l) as {r*, C Xfc = tfct^}. 
By fixing the transmit covariance matrices as {C Xfc = tfctj^}, 
(PI) is further optimized over the pseudo-covariance matrices 
{C Xfc } in this subsection. By replacing the first term in the rate 
expression ([8) by a k r* , the problem for pseudo-covariance 
matrix optimization is formulated as 



(PI. 3): max 

fl,{C x .} 



R 



s.t. a k r 



-log- 

2 S 1-C S - 2 |C.. 



> a k R, Vfc, 



tfct 



(t fc t£r 



>- 0, Vfc, 



(11) 



where C Vk and C Sk are fixed covariances given the previously 
optimized transmit covariance matrices {C Xfc = tfct^}; C Vk 
and C Sk are the pseudo-covariances, which are related to 
transmit pseudo-covariance matrices {C Xfc } via Q and ((71, 
respectively. By treating R as a slack variable and discarding 
irrelevant terms in (PI. 3), the problem can be re- written as 



a minimum-weighted-rate maximization (MinWR-Max) prob- 
lem as follows. 



(PI. 4): max. 
l c * fc } 

s.t. 



1 , 1 

mm . log — 

fc=l,- ,K 2(1/, 1 



C y k 2 \ C Vk\ 2 

C7 2 \C Sk \ 2 



tfctf 



(tfctf)* 



y o, Vfc. 



(12) 



(PI. 4) is a matrix optimization problem that deals with the 
additional rate term in ([8]) due to the use of improper Gaussian 
signaling. Denote the optimal value of (PI. 4) as r*. If r* > 0, 
then a strict rate improvement corresponding to rate-profile a 
over the optimal proper Gaussian signaling is achieved. The 
following result will be used for solving (PI. 4). 

Lemma 1. The positive semidefinite constraint in \\2\ is 
satisfied if and only if 

C Xt = Z k t k tl, Vfc, (13) 

where Z k is a complex scalar variable satisfying \Z k \ < 
||t/j|| 2 , and tfc = tfc/||tfc|| denotes the normalized transmit 
beamforming vector for user fc with proper Gaussian signal- 
ing. 

Proof: Please refer to Appendix [A] ■ 
It is noted that C Xfc in ( fT3| ) is a rank-1 symmetric matrix, 
and thus Lemma [T] establishes the optimality of rank-1 pseudo- 
covariance matrices in the if -user MISO-IC, if the optimal 
rank-1 transmit covariance matrices are applied. Furthermore, 
it follows from ( fT3| l that each rank-1 pseudo-covariance matrix 
is parameterized by one single complex scalar Z k . As a result, 
the problem dimension for pseudo-covariance optimization can 
be significantly reduced from KM 2 in the original matrix 
problem (PI. 4) to K by applying Lemma [T] as will be shown 
next. 

By substituting ( fT3j ) into |5]) and |7ji, we have 

K 

c Vk = ^(hk&fzj, c Sk = >>, ; l;-z,. vfc. (14) 



Define the following if -dimensional complex-valued vectors: 
z = [Zi ■ ■ ■ Z K ] T , 
mfc = Cy k x [(hfeiti) 2 

w fe = C s k [• ■ ■ ( h k(k 
Then we have 



• • • (hfciftx) 2 ] , 
-i)tfc-i) 2 (h fc ( fe+1 )t fc+1 ) 2 



c y k \ c yh\ 

cz 2 \c Sk \ 

Sk I °K | 



ifz| 2 , 



Iw^zl 2 . 



(15) 
(16) 



Therefore, (PI. 4) can be reformulated as 

1 , 1 - |m 



(PI. 5): max . 

zec K 

s.t. 



mm . los 

fc=i,- ,K 2a k 



H z\ 2 



1 - |wf zP 



< tt 



, Vfc, 



(17) 



where e k is the fcth column of a K x K identity matrix. 
Note that the constraint in ( fT7] i corresponds to the condition 
\Z k \ < || tfc || 2 given in Lemma [T] Before solving (PI. 5), we 
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first give an intuitive discussion on when it is possible for 
(PI. 5) to have a strictly positive objective value, or in other 
words, achieve a strict rate improvement over the optimal 
proper Gaussian signaling by further optimizing the pseudo- 
covariance matrices. 

Denote the second term of the rate expression in <|8j as 
ARk, which is the additional rate gain due to the use of 
improper Gaussian signaling. It is easy to verify that with the 



optimal proper Gaussian signaling obtained in Section III-A 



with C,, 



0, Vfc, we have ARk = 0, Vfc. In contrast, with 



the rank-1 transmit co variance and pseudo-covariance matrices 
given above, a simple upper bound on ARk for (PI. 5) can be 
obtained as 



AR k 



1 i-K'z 2 

2 l-|w£z| 2 

i, i 

- 1Og l 1 H 12 

2 1 — |w" zp 

<-ilog(l-||w fe || 2 ||z|| 2 ) 



(18) 
(19) 
(20) 



where in order for the upper bound in ( |20| > to be meaningful, 
||wfc|| 2 is assumed to be small enough so that ||wfc|| 2 ||z|| 2 < 1. 
Note that ||w fe || 2 = C~ L 2 J2j^k l^-kftjl 4 , which reflects the 
interference level experienced at receiver fc given the beam- 
forming vectors {tj.}. As ||wj.|| 2 — > 0, the upper bound in 
( |20| i approaches to 0, which leads to the following remark. 

Remark 1. Improper Gaussian signaling is more advan- 
tageous than proper Gaussian signaling in MISO-IC only 
when there is non-negligible interference among the users. 
For example, with zero-forcing (ZF) transmit beamforming 
vectors {tf F }, i.e., h^tf^ = 0> ^3 h we have ||w fe || 2 = 
and thus ARk = 0, Vfc. As a result, no rate improvement 
can be achieved by further optimizing the pseudo-covariance 
matrices. This is as expected since ZF transmit beamforming 
essentially results in K independent point-to-point MISO 
channels, where proper Gaussian signaling is known to be 
optimal ^38^ . 

For the general K-usei MISO-IC, ZF transmit beamforming 
is feasible only when the number of transmitting antennas 
at each transmitter is no smaller than the total number of 
users, i.e., M > K. For the general case where the users are 
interfered by each other with non-ZF transmit beamforming, 
i.e., Wfe ^ 0, there is then a potential opportunity to improve 
the achievable rates over the optimal proper Gaussian signaling 
by solving the pseudo-covariance optimization problem in 
(PI. 5), as we pursue next. 

(PI. 5) is a problem of maximizing the minimum of K 
weighted log-fractions of quadratic functions over z, for which 
the well-known SDR technique can b e applied to find an 
efficient approximated solution |32 33j. Let Mj. = m^m^ 
and Wfc = WfcW^, Vfc, and with the identity x H Ax = 
Tr(Axx ff ), the SDR problem of (PI. 5) is given by 



(P1.5-SDR): max. min. — log- — Tr (^ Ifc ^) 
zbo k=i,-,K 2a h 6 l-Tr(W fe Z) 



where Efe = e^e^ . It is easy to see that (PI. 5) is equivalent to 
(P1.5-SDR) with the additional constraint rank(Z) = 1, under 
which Z can be written as Z = zz H . Therefore, the optimal 
value of (P1.5-SDR), denoted as r s d r , provides an upper bound 
on that of (PI. 5). Since Z = is always feasible for (PI. 5- 
SDR), which corresponds to an objective value of 0, we have 
Tsdr > 0. In other words, the proposed pseudo-covariance 
optimization is guaranteed to improve the achievable rates 
over the optimal covariance (beamforming) optimization in 
Section IiTTa! 

Theorem 1. For any matrix Z that is feasible to (P1.5-SDR), 
the following inequalities hold: 

1 - Tr(W fc Z) > C~V > 0, Vfc, (22) 
1 - Tr(M fc Z) > C^ct 4 > 0, Vfc. (23) 

Proof: Please refer to Appendix B] ■ 
With the inequalities in p2} , ( |23 1 and r s( j r > 0, it then 
follows that (P1.5-SDR) is a quasi-convex problem p7] . For 
any given r > 0, we consider the following problem. 

(P1.6): min. Tt^Z) 



s.t. 1 - Tr(M fc Z) > e 
Tr(E fc Z) < ||t fe || 4 



2QfcT (l-Tr(W fc Z)), Vfc, 
fc = 2,--- .K. 



(PI. 6) is a semidefinite programming (SDP) problem, which 
minimizes the left hand side (LHS) of ( pTj l corresponding to 
fc = 1, subject to a target objective value r for (P1.5-SDR). 
Denote the optimal value of (P1.6) as /(r). If /(r) < ||ti|| 4 , 
then the optimal value of (P1.5-SDR) satisfies r s d r > r; 
otherwise, r s d r < r. Therefore, (P1.5-SDR) can be optimally 
solved by solving the SDP problem (PI. 6) with different values 
of r, and applying a bisection search over r. 

Denote the solution to (P1.5-SDR) by Z*. If rank(Z*) = 1, 
i.e., Z* = zz H , then z is the optimal solution to (PI. 5). In this 
case, our proposed SDR is tight; otherwise, if rank(Z*) > 1, 
we then apply the following Gaussian randomization pro- 
cedure customized to our problem to find an approximated 
solution to (P1.5) (32). 

Algorithm 1 Gaussian Randomization Method for (PI. 5) 
Input: The solution Z* to (P1.5-SDR), and the number of 
randomization trials L. 

1: for / = 1, • • • ,L do 

2: Generate ~ C7V(0, Z*), and construct a feasible 
point z; to (PI. 5) as follows: 



( lit II 2 ~i 

[zi]k = Kfc[&] fc ,with K k = minjl, . j, 



end for 

determine I* 



arg max . 



1=1,— ,l fc=i,— ,K 



Vfc. 



zfM fc z, 



z"W fc z, 



s.t. Tr(EfcZ) < lit* 



Vfc, 



(21) 



Output: z = z;* as an approximated solution for (PI. 5). 

To summarize, our proposed separate covariance and 
pseudo-covariance optimization algorithm for (PI) is given in 
Algorithm [2] 
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Algorithm 2 Separate Covariance and Pseudo-Covariance 
Optimization for (PI) 

Solve (Pl.l), denote the solution as {r* , C* fc = t k tf }. 
Solve (P1.5-SDR), denote the solution as Z*. 
If rank(Z*) = 1, then its principal component z with 



Z* = zz is the optimal solution to (PI. 5); otherwise, 
find an approximated solution z using Algorithm [T] 
Obtain the pseudo-covariance matrix solution {C* fc } us- 
ing ([13), and the maximum sum-rate R* using ([8]). 



Remark 2. In the special case of K — 2, (PI. 6) is a complex- 
valued SDP problem with three affine constraints. It is known 
that if such a problem is feasible, there is always a rank- 
1 optimal solution j |59| /. Therefore, for the two-user MISO- 
IC with rank-1 transmit covariance matrices, the SDR-based 
solution will give the optimal pseudo-covariance matrices of 
rank-1. Note that in our previous work \31^ , a SOCP-based 
algorithm has been proposed, which is able to find the optimal 
pseudo-covariances for the two-user SISO-IC. However, the 
SOCP-based algorithm is difficult to be extended to the case 
of K > 2 and/or MISO-IC with M > 1. 

IV. Improper Gaussian Signaling for MISO-BC 
with Linear Precoding 

The improper Gaussian signaling scheme discussed in the 
preceding sections can also be applied to the Gaussian MISO- 
BC under the assumption of linear precoding employed at the 
base station (BS) transmitter. Consider a Gaussian MISO-BC 
where a transmitter with M > 1 transmit antennas sends 
independent information to K single-antenna receivers. The 
signal received by the fcth user can be written as 



Vk = hfcx + v k 

= h fc x fe + ^2 h k x j + v k, k = l, 



(24) 



where h& G C denotes the channel vector from the 
transmitter to user fc; v k ~ CJ\f(0,a 2 ) represents the CSCG 
noise; and x = Y2k=i Xfc * s ^ e transmitted signal vector 
from the transmitter with x& denoting the transmitted signal 
intended for user fc. 

It is known that the capacity of the MISO-BC can be 
achieved by employing the "dirty paper coding (DPC)" tech- 
nique at the transmitter [40]. DPC is a nonlinear precoding 
technique by which the information for different users is 
encoded in a sequential manner, so that the interference caused 
by earlier encoded users can be completely removed at the 
transmitter with the non-causal information of the earlier 
encoded users' data streams. Since the nonlinear capacity- 
achieving DPC scheme is difficult to implement in practical 
systems, a great deal of research has focused on simpler linear 
precoding schemes at the transmitter. In this case, all the 
inter-user interferences are treated as additive Gaussian noise, 
thus resembling a MISO-IC (but with the transmitter's sum- 
power constraint as opposed to transmitters' individual power 
constraints in the general MISO-IC defined in Section 



Denote by t k the transmit beamforming vector for user k in 
a K-usei MISO-BC with linear precoding, we then have 



Xfc = tfedfc, Vfc, 



(25) 



where d k is a CSCG random variable representing the in- 
formation signal of user k, i.e., d k ~ CAT(0, 1). The signal 
vector Xfc in ( |2"5j ) is a proper Gaussian RV with covariance 
and pseudo-covariance matrices respectively given by 

C Xfc = t fc tf , C Xfc = 0, Vfc. (26) 

However, it is not immediately clear from |24) whether the 
restriction of proper Gaussian signaling, or zero pseudo- 
covariance matrices of the transmitted signals, will incur any 
rate loss under the assumption of linear precoding in the 
MSIO-BC. Therefore, we consider the more general improper 
Gaussian signaling similar to the MISO-IC, where the pseudo- 
covariance matrices {C Xfc } are further optimized. Following 
the similar derivation as for the MISO-IC in ([8j, the achiev- 
able rate for user fc in a MISO-BC with improper Gaussian 
signaling can be expressed as 



Rk 



BC 



log 1 



E 



j7tfc |hktj| 2 



+ W"^ 21 -^ 2 (27) 

where C yk , C Vk , C Sk , and C Sk are defined similarly as in 
(|4|-(|7|. Note that the rate expression ( f27] ) bears the same 
structure as ([H) for the MISO-IC. For a given rate-profile 
a.= (ai ■ ■ ■ ax), the problem to find a Pareto-optimal rate- 
tuple for the MISO-BC can thus be formulated as 

(Pl-BC): max. R 



s.t. R k 3C > a k R, Vfc, 



H 



(tfct 



x/, 
H\* 
k ) 



K 

Ei 

k=l 



h 0, Vfc, (28) 



(29) 



where |28) is the necessary and sufficient condition for t^tf 
and C Xfc to be a valid pair of covariance and pseudo- 
covariance matrices, and (|29) denotes the sum-power con- 
straint at the transmitter. Note that from the optimization 
perspective, (Pl-BC) is identical to (PI), and hence the sepa- 
rate covariance and pseudo-covariance optimization algorithm 
presented in the preceding section for the MISO-IC can be 
directly applied to solve (Pl-BC). 

V. Numerical Results 

In this section, we evaluate the performance of the proposed 
algorithm by numerical examples. All the transmitters are 
assumed to have the same power constraint P, i.e., Pf. = 
P, Vfc. The average SNR is defined as P/a 2 . For the Gaussian 
randomization in Algorithm [T] L = 1000 is used. 
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TABLE I: Mean and standard deviation (std) of the approximation 
ratio upper bound t^/t. 





K 


2 


3 


4 


5 


6 


M = 1 


mean 


1.0 


1.032 


1.138 


1.267 


1.391 


std 





0.092 


0.245 


0.350 


0.441 


M = 2 


mean 


1.0 


1.012 


1.162 


1.401 


1.640 


std 





0.068 


0.388 


0.621 


0.691 



— A — Improper, SDR-based separate opt. 
— 9 — Optimal proper 

























1 1.5 
R, (bits/channel use) 



Fig. 1: Achievable rate region for the two-user MISO-IC with 
M = 2, SNR = 10 dB, and channel realization H (1) . 



A. Approximation Ratio for SDR 

First, we evaluate the quality of the SDR-based approx- 
imated solution for the pseudo-covariance optimization pro- 
posed in Section III-B Denote r* and r sc j r as the optimal 



objective values of (PI. 5) and its relaxation (P1.5-SDR), 
respectively. Also denote f as the objective value of (PI. 5) 
with the approximated solution obtained by Algorithm [T] Then 
we have 

f < t* < T sdr , or 1 < t*/t < T sdl -/f , 

where r*/f is the approximation ratio. Since in general the 
optimal value t* is difficult to be found, the upper bound 
Tsdr/'r of the approximation ratio can be used to evaluate 
the quality of the SDR-based approximated solution (33}. If 
7" s dr I 't = 1, then the obtained SDR-based solution is in fact 
optimal. With the rate-profile in (PI) given as a = 1/K1, 
where 1 is an all-one vector, Table [I] summarizes the mean and 
standard deviation of t^/t at SNR = 10 dB with different 
pairs of values for M and K over 1000 random channel 
realizations, each with the channel coefficients drawn from 
the i.i.d. CSCG random variables with zero-mean and unit- 
variance. It is observed that for all the setups considered, 
the mean values of the approximation ratio upper bounds are 
between 1 and 1.64. In particular, for K = 2, r s( j r /f = 1 is 
guaranteed, which verifies the optimality of the SDR-based 
pseudo-covariance optimization for the two-user MISO-IC, as 
discussed in Remark [2] 



- Improper, SDR-based separate opt 

- Optimal proper 




1 1.5 
R, (bits/channel use) 



Fig. 2: Achievable rate region for the two-user MISO-IC with 
M = 2, SNR = 10 dB, and channel realization H (2) . 



column of Table [IT] The proposed improper Gaussian signaling 
with SDR-based separate covariance and pseudo-covariance 
optimization is compared against the optimal proper Gaussian 
signaling obtained by solving (PI. 1). It is observed that for this 
channel setup, the achievable rate region has been significantly 
enlarged by the proposed pseudo-covariance optimization for 
improper Gaussian signaling. 

Next, we consider a different channel realization H^ 2 ', 
which differs from H' 1 ) only in the phases of the second 
elements in and hi2, as shown in the right column of 
Table |H| It is observed from Fig. [2] that for this new channel 
setup, there is no notable rate gain by the proposed improper 
Gaussian signaling over the optimal proper signaling, which is 
in sharp contrast to the result in Fig. [T] This can be explained 
by comparing the residue interference levels after applying 
the optimal proper Gaussian signaling in the two cases. It 
is worth noting that for the two-user MISO-IC, the optimal 
proper Gaussian signaling at each transmitter is known to be 
a linear combination of the ZF beamforming and the maximal 
ratio (MR) beamforming [6|. Applying this result to the two- 
user MISO-IC in our context, it then follows that one user 
will generate less amount of interference to the other user if 
its direct channel and interfering channel vectors are more 
close to be orthogonal]^] Let B\ denote the angle between 
the direct and interfering channel vectors for transmitter 1, 
i.e., cos^i = |h2ih^|/(||h2i||||hn||), and define #2 for 
transmitter 2 similarly. As shown in Table |II} since 9\ and O2 
are smaller in the case of than that in the case of K^ 2 \ 
more substantial interference is resulted even after applying 
the optimal proper transmit beamforming. As a result, further 
optimization over the pseudo-covariance matrices provides 
more notable rate gains in the case of than HS 2 ', as 

discussed in Remark Q] 



B. Rate Region Comparison 

In Fig. [T] the achievable rate region for an example two- 
user MISO-IC with M = 2 is plotted for SNR=10 dB, with 
the channel realization (denoted as H' 1 )) given in the left 



2 Consider one extreme case when the direct and interfering channel vectors 
are orthogonal, then the MR and ZF beamforming vectors coincide, and thus 
no interference is generated to the other user; however, for the other extreme 
case when the direct and interfering channel vectors are linearly dependent, 
a nonzero interference is always resulted. 
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TABLE II: Channel realizations for Figs. [TJ and [2] 





channel realization H^ 1 ' 


channel realization H 1 - 2 ' 




[0.3676e- 1 - yu;iYl 0.4993e 1 - bUYbl ] 


[0.3676e- i - YU3Yl 0.4993e ibu ' bl ] 


h 2 i 


[0.2526e- 1SMaYl 0.3270e 1&SB41 ] 


[0.2526e- i!,ayYs 0.3270e u a81Ul ] 


h 2 2 


[0.4694e- u - 19i! " 0.5682e u - btlU21 ] 


[0.4694e- u - iylbI 0.5682e u ' bilU21 ] 


hi2 


[0.2885e- U24541 0.3656e u 47lul ] 


[0.2885e- u - 24548 0.3656e lsli7ai ] 


cos 9i 


0.9961 


0.6601 


cos 8 2 


0.9997 


0.7793 




Fig. 3: Average max-min rate for the two-user SISO-IC. 



— ©— Improper, SDR-based separate opt. 
— * — Optimal proper 





































15 

SNR (dB) 



Fig. 4: Average max-min rate for the three-user MISO-IC with 
M = 2. 



C. Max-Min Rate Comparison 

The rate-profile technique used in characterizing the Pareto 
boundary of the achievable rate region can be directly applied 
for maximizing the minimum (max-min) rate of all the users 
without time sharing. Specifically, the max-min problem for 
the X-user MISO-IC is equivalent to solving (PI) by using the 
rate-profile ol — 1/K1. In this subsection, the average max- 
min rate achievable over 1000 random channel realizations by 
the proposed improper Gaussian signaling is compared with 
that by the optimal proper Gaussian signaling obtained by 
solving (Pl.l). 

We first consider the special case of the two-user SISO-IC, 
where two other existing improper Gaussian signaling designs 



based separate covariance and pseudo-covariance optimiza- 
tion, are compared with the proposed SDR-based separate 
covariance and pseudo-covariance optimization. It is observed 
in Fig.[3]that improper Gaussian signaling provides significant 
gains over the optimal proper Gaussian signaling. At high 
SNR, the max-min rate by proper Gaussian signaling saturates, 
since the total number of data streams transmitted, which is 
2, exceeds the total number of DoF of the two-user SISO-IC, 
which is known to be 1. In contrast, a linear increase of the 
max-min rate with respect to the logarithm of SNR is observed 
by improper Gaussian signaling. It is also observed that the 
two separate covariance and pseudo-covariance optimization 
algorithms, SDR-based or SOCP-based, provide the same 
max-min rate in the case of K = 2. This is as expected 
since both algorithms achieve the global optimality for the 
covariance and pseudo-covariance optimization subproblems 
in the two-user SISO-IC case. Moreover, it is observed in 
Fig. [3] that the SDR-based joint covariance and pseudo- 
covariance optimization algorithm in |31| achieves additional 
rate gains over the SDR/SOCP-based separate optimization. 
However, the extension of such a joint optimization to the 
general X-user MISO-IC remains unknown. 

To show the max-min rate performance with improper Gaus- 
sian signaling when there are multiple transmitting antennas, 
we consider a three-user MISO-IC with M = 2. As shown 
in Fig. [4] a significant rate improvement is observed by the 
proposed improper Gaussian signaling over the optimal proper 
Gaussian signaling. 



D. Rate Region of MISO-BC 

At last, MISO-BC with linear precoding as discussed in 
Section IV is considered. In Fig. [5] the achievable rate regions 
for an example two-user MISO-BC with M = 2 are plotted 
for SNR = 10 dB. The channel vectors of the two users are 
respectively given by hi = [1.1741e 1 0030i 0.8064e 2 86421 ], 
h 2 = [1.8116e 2 0647i 0.9209e- 2 ' 41671 ]. Under the assump- 



proposed in |31|, one with SDR-based joint covariance and 
pseudo-covariance optimization and the other with SOCP- 



tion of linear precoding, the proposed improper Gaussian 
signaling is compared with the optimal proper Gaussian sig- 
naling with linear beamforming. As a benchmark, the capacity 
region achieved by DPC is also included. It is observed that 
without time sharing (TS) of the achievable users' rates by 
a convex-hull operation, the achievable rate region for the 
MISO-BC is significantly enlarged by employing improper 
Gaussian signaling over the optimal proper Gaussian signaling. 
However, such a performance gain no longer exists if TS is 
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— - Dirty Paper Coding 

— A — Improper, SDR-based separate opt. 

— © — Optimal proper 

Time sharing 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 
R 1 (bits/channel use) 

Fig. 5: Achievable rate region for the two-user MISO-BC 
employing linear precoding with M = 2, and SNR = 10 dB. 



performed^] 



VI. Conclusion 



This paper studies the transmit optimization problem for 
the if -user MISO-IC with the interference treated as Gaus- 
sian noise, when improper or circularly asymmetric complex 
Gaussian signaling is applied. By exploiting the separable 
achievable rate structure with improper Gaussian signaling, a 
separate transmit covariance and pseudo-covariance optimiza- 
tion algorithm is proposed. For the pseudo-covariance opti- 
mization, we show the optimality of rank-1 pseudo-covariance 
matrices, given the optimal rank-1 transmit covariance matri- 
ces obtained by existing methods. Moreover, we show that 
each rank-1 pseudo-covariance matrix is parameterized by 
one single complex scalar and thereby greatly reduce the 
dimension for searching. We then apply the SDR technique 
to find an efficient approximated solution for the pseudo- 
covariance optimization. The proposed improper Gaussian 
signaling is also extended to the case of MISO-BC with 
linear precoding. Finally, numerical results are provided to 
demonstrate the achievable rate gains over the conventional 
proper Gaussian signaling in multiuser MIMO systems that 
are modeled by the MISO-IC. 



We now re-express the positive semidefinite constraint in 
( fT~2| ) using the above three conditions in Lemma [2] First, it is 
clear that (a) is guaranteed in \V2) . Next, to use the condition 
in (b), we express the singular value decomposition (SVD) of 



the rank-1 covariance matrix CJ as 



C i fc _ tkt k — [tfc Ufc] 



tfc 



z k 



where t k = t k /\\t k \[, and U fc e cAfx(M-i) satisfies 
Uf U fe = Im-i. and tf U fe = 0. Then the Moore-Penrose 
pseudo-inverse of C~^ k can be obtained as 

(c:jt = nt fc ir 2 t fc tf. 

Applying the condition in (b) of Lemma [2] to ( fl2] > yields 
(I - C* (C* fc )t) C Xfc = (I - t fc tf )C X , = 



7*tfcv£ 



H 



(30) 



where j k is a nonnegative real number denoting the singular 
value of C Xfc , and v k G C M is a unit-norm complex vector. 
Since C Xfc is a pseudo-covariance matrix, which must be 
symmetric, we have 



By expressing C Xk C^ using two alternative forms, we have 

tfc = e""vl, or v fe = e 
where 9 k & [0, 2ir). By substituting v k into ( (30] ), we have 



i8 "v* k , or v fe = e ie "t* k , (31) 



7fee 



e tfct 



fe L fc 



(32) 



where we have defined the complex variable Z k = j k e~ l6k . 
Furthermore, the condition in (c) of Lemma [2] implies that for 
the constraint in ( fT2] > to be satisfied, we need to have 



c* fe (c x jtc x ,ho 



(iit fc || 2 -|z fe | 2 nt fe ir 2 )t*t^o 

\Z k \ < ||tfe|| 2 . 
This completes the proof of Lemma [T] 



(33) 



Appendix A 
Proof of LemmaQ] 

The following result will be used for proving Lemma [T] 

Lemma 2. [37] If X is Hermitian and is partitioned as 

X = —if _ , then X if and only if the following 

three conditions are satisfied: 

(a) A y 0; (6) (I - AA*)B = 0; (c) C - B^A^ h 0, 
where (-)t represents the Moore-Penrose pseudo-inverse. 

3 Based on our simulations over many different channel and SNR setups 
(not shown here due to the space limitation), it is found that the rate gain 
by the improper signaling over the optimal proper signaling diminishes after 
the TS is applied for the MISO-BC with linear precoding, for which some 
further investigation is needed to reveal its fundamental reason. 



Appendix B 
Proof of Theorem[TJ 

Since the inequalities in ( p2) and ( |23] l of Theorem [T] can be 
proved similarly, for brevity, we only show the proof of p2| 
in this appendix. First, if = w k wjf = 0, i.e., w k is a 
zero-vector, then d22| is satisfied since 



a 2 ?>a 4 



(34) 



is true. Therefore, in the following, we assume without loss of 
generality that at least one of the elements in w k is non-zero. 
Then we consider the following optimization problem. 



(P-Al): max. Tr(W fc Z) = Tr(w fc wf Z) 
s.t. Tr(E fc Z) < lltJI 4 , Vfc. 



(35) 



to 



In order to show |22) in Theorem [T] it is sufficient to prove 
that the optimal objective value of (P-Al) is no greater than 
1— C,T 2 <7 4 . With ([34]i, this is clearly true if the optimal solution 



Z* to (P-Al) is a zero matrix. Therefore, in the following, we 
assume Z* 7^ 0. First, we give the following lemma. 

Lemma 3. There exists a rank-1 optimal solution to (P-Al). 

Proof: Lemma [5] can be proved by applying the more 
general result given in Lemma I of | |4T| , where Ej. can be any 
positive semidefinite matrix. For the problem of our interest 
where in ( |35| ) is a matrix with all zero-elements except 
that the (k, fc)-th entry equals to one, we provide an alternative 
proof for the same result based on the Karush-Kuhn-Tucker 
(KKT) optimality conditions to provide more insight. Without 
loss of generality, by rearranging terms if necessary, we can 
partition Wfc 6 C K as 



Wfc 





(36) 



where Wfc £ C A consists of all non-zero elements in 
where 1 < K' < K. Accordingly, we partition Z € 



-<KxK 



as 



Z Z12 

Z21 Z22 



(37) 



where the dimensions of Z, Z12, Z21 and Z22 are given by 

K' x K', K' x (K ~ K'), (K - K') x K' , and (K - K') x 
(K — K'), respectively. Then we have 



w fe Zwfc = w fe Zwfc. 



(38) 



In other words, the matrices Z12, Z21 and Z22 are not related 
to the objective value of (P-Al). Furthermore, since Efc in ( f3~5| ) 
is diagonal, without loss of optimality, Zi2,Z2i and Z22 can 
all be set to zero matrices. Then (P-Al) can be recast as 



(P-A2): max. Tr(w fe wfZ) 
z 

s.t. Tr(E fe Z) < ||t fc | 

Z h 0, 



(39) 
(40) 



where E& e 



is a matrix with all zero-elements except 



that the (fc, k)-th entry equals to one. Next, we show that the 
optimal solution to (P-A2) is rank-1 by investigating its KKT 
conditions. Specifically, let Z* be the optimal solution to (P- 
A2), {X* k > 0} and Y* y be the optimal dual solutions 
associated with the constraints ( |39"1 > and ( |40| ), respectively. 
According to the KKT conditions of (P-A2), one can verify 
that the following conditions must be satisfied: 



Y* = A - w fc wf h 0, 



(41) 
(42) 



Y*Z* = 

where A = £f=i A£E fc = Diag{A}, • • • ,\* K ,}. From (|4T> : 
we have 



A*>0, * = !,...,#', 



(43) 



which can be shown as follows by contradiction. Suppose that 
there exists an index I e {!,■■■ ,K'} with = 0. Then we 



define a vector v = 5eu where 5 > and e; is the Ith column 
of the identity matrix \k>- We thus have 

v H Y*v = -5 2 \e? w fc | 2 = -(^[w^l 2 < 0, (44) 

where the last inequality follows since all the elements in 
Wfe are non-zero. ( p4| ) clearly contradicts with the positive 
semidefinite constraint of Y* in ( |4T| , and hence ( |4"3j ) must be 
true. Then we have A = As As, with As being invertible. 
Thus, we have 

rank(Y*) = rank ^A^ \L K , — A^^wj-wf A _ 5]A5 

= rank (Ik 1 — A _3 WfeW ^A _ s 
>K' — 1. (45) 
It follows from ( |42| i and (j45j that Z* must be rank-1 since 
< rank(Z*) < Nullity (Y*) = K' - rank(Y*) < 1. (46) 



Given f37] i and the zero matrices for Z12, Z21 and Z22, the 
solution Z* to (Pl-Al) must be of rank-1 since 



rank(Z*) = rank 



Z* 




rank(Z*) = 1. (47) 



This completes the proof of Lemma [3] 

With Lemma [3] (P-Al) can be re-expressed as 

(P-A3): max. z H W k z 



s.t. |e?z| 2 < ||ti,|| 4 , 



Vfc. 



(48) 



Using ( fT4l > and ( p~6] >, we have 
z H W k z 



a- 2 IE( h ^. 



■-1. ! / 



<C- 2 (^|h fei t,| 2 |^ 

(^Ih^fllt.H 2 



< c; k 2 



3^k 
-2 I / \ " 



<C-[( 



= 1 - 



(49) 

(50) 

(51) 
(52) 



where ( |49| > is due to the triangle inequality; ( |50j > is due to the 
constraint in (j48j, which is equivalent to \Z k \ < ||tfc|| 2 ,Vfc; 
and ( |52"] > is true since C Sk = X^j^fc l^fcj^il 2 + a2 - The above 
result shows that the optimal objective value of (P-A3), and 
hence that of (P-Al), is no larger than 1 — C^cr 4 , as desired. 
This completes the proof of Theorem [T] 
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